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Abstract 

Recently presented explicit formulae for asymptotic expansions of Feyn- 
man diagrams in the Sudakov limit [|l| are applied to typical two-loop 
diagrams. For a diagram with one non-zero mass these formulae provide 
an algorithm for analytical calculation of all powers and logarithms, i.e. 
coefficients in the corresponding expansion (Q 2 )~ 2 J2 n ,j=o c njt~ n In- 7 1, 
with t = Q 2 /m 2 and j < 4. Results for the coefficients at several first 
powers are presented. For a diagram with two non-zero masses, re- 
sults for all the logarithms and the leading power, i.e. the coefficients 
c n j for n = and j = 4,3,2,1,0 are obtained. A typical feature of 
these explicit formulae (written through a sum over a specific family 
of subgraphs of a given graph, similar to asymptotic expansions for 
off-shell limits of momenta and masses) is an interplay between ultra- 
violet, collinear and infrared divergences which represent themselves 
as poles in the parameter e = (4 — d)/2 of dimensional regularization. 
In particular, in the case of the second diagram, that is free from the 
divergences, individual terms of the asymptotic expansion involve all 
the three kinds of divergences resulting in poles, up to 1/e 4 , which are 
successfully canceled in the sum. 
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1 Introduction 



The simplest explicit formulae for asymptotic expansions of off-shell Feynman dia- 
grams in various limits of momenta and masses^ 0. [5], |3j (see || for a brief review) 
have been recently generalized jl], ^| for two on-shell limits: the limit of the large 
momenta on the mass shell with the large mass and the Sudakov limit, with the large 
momenta on the massless mass shell. 

To derive these formulae a method based on the construction of a remainder of 
the expansion and using then diagrammatic Zimmermann identities [[J and applied in 
refs. || for operator product expansions within momentum subtractions and later for 
diagrammatic and operator expansions within dimensional regularization and renor- 
malization |4j] has been used. These general prescriptions were illustrated in ref. [1] 
through one-loop examples (the triangle diagram of Fig. la in the case of the Sudakov 
limit) and, in ref. |J, illustrated and applied for calculation of a typical two- loop di- 
agram in the case of the first of two limits. 

The purpose of this paper is to present similar results in the latter case, i.e. for the 
Sudakov limit. We shall consider two typical vertex two- loop diagrams. The former 
is the Feynman integral F 1 (p 1 ,p 2 , m) corresponding to Fig. 2a with mi = . . . = m 5 = 
0, m 6 = m, and the latter is F 2 (p 1 ,p 2 , tti) with mi — . . . — m± — 0, m 5 = m 6 = m. 
The limit under consideration is q 2 = (pi — p 2 ) 2 = —Q 2 — ► — oo, with p\ = p\ = 0. 
In the case of Fx, the general prescription will be used to obtain an algorithm for 
analytical calculation of arbitrary coefficients c£9 in the corresponding expansion 

„2 _ QO oo 4 

F l (pi, P2 ,m) ~ (Q 2 )- 2 J2J:^(m 2 /Q 2 r\^(Q 2 /m 2 ). (1) 

71=0 j = 

Results for the coefficients at several first powers will be presented. For F 2 , the 
coefficients for all the logarithms and the leading power, i.e. c^- for n = and 
j = 4, 3, 2, 1, will be analytically evaluated. 

As in the case of the off-shell explicit formulae ^ |3|, || the prescriptions of ref. [IJ 
are written through a sum over a specific family of subgraphs of a given graph. 
A typical feature of the off-shell formulae is an interplay between ultraviolet and 
infrared divergences which appear in individual terms of that sum but are mutually 
canceled, provided the initial Feynman integral is finite. It turns out that, for the 
Sudakov limit, one meets a more general interplay between ultraviolet, collinear and 
infrared divergences which, at first sight, seem to be of very different nature. In 
particular, in the case of F 2 , which does not have divergences from the very beginning, 
individual terms of the corresponding asymptotic expansion involve all three kinds of 
divergences resulting in poles up to 1/e 4 , with e = (4 — d)/2 parameter of dimensional 
regularization ||. However these poles are successfully canceled in the sum and one 
obtains expansion (Q) in the limit e — > 0. 

2 In the large mass limit, one can also apply the same 'off-shell' formulae for any values of the 
external momenta. 
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In the next section we shall briefly recall the general prescriptions and one-loop 
example of ref. 0. In the following sections these prescriptions will be applied to the 
Feynman integrals F\ and F2. 



2 General prescriptions and one-loop example 

The asymptotic expansion of a Feynman integral Fr(pi,p 2 ,m), corresponding to a 
vertex graph T, in the Sudakov limit takes the following explicit form Q: 

q 2 — > —00 

F T (p 1 ,p2,rn,e) «~ 2^ M^Fr(pi,P2,m, e) . (2) 

7 

Here the sum runs over subgraphs 7 of T for which at least one of the following 
conditions holds: 

(i) In 7 there is a path between the end-points 1 and 3. (Three end-points of the 
diagram are numerated according to the following order: pi,P2,Q = Pi — P2-) The 
graph 7 obtained from 7 by identifying the vertices 1 and 3 is 1PI. 

(ii) Similar condition with 1 <-> 2. 

The pre-subtraction operator M.^ is defined as a product WjM. lj of operators of 
Taylor expansion acting on 1PI components and cut lines of the subgraph 7. Suppose 
that the above condition (i) holds and (m) does not hold. Let jj be a 1PI component 
of 7 and let p\ + k be one of its external momenta, where k is a linear combination 
of the loop momenta. (We imply that the loop momenta are chosen in such a way 
that pi flows through all jj and the corresponding cut lines). Let now q. be other 
independent external momenta of jj. Then the operator M. for this component is 
defined as % ! -((p 1 k)/(pip 2 ))p2 > q.,m- > where rrij are the masses of jj. In other words, it 
is the operator of Taylor expansion in q. and rn.j at the origin and in k at the point 

k = ^^\ p2 (which depends on k itself). 

For the cut lines the same prescription is adopted. If pi + k is the momentum 

of the cut line, then the corresponding operator acts as T K K ^_^ +2pik 

conditions (i) and (ii) hold the corresponding operator performs Taylor expansion in 
the mass and the external momenta of subgraphs (apart from p\ and P2). 
For example, in the case of the one-loop triangle diagram of Fig. la, 

F r ( Pl ,p 2 , m,e)=J {k ,_ 2vik) ^_ k 2V 2kW-m-) ' (3) 

this family of subgraphs is shown in Fig. 1. Besides the graph itself (a), these are two 
subgraphs consisting, respectively, of the left and the right line (b and c). The sub- 
traction operator for (a) expands the propagator with the mass m in Taylor series in 



. If both 

1 



2 



m, and the subtraction operators A4 7 for (b) and (c) act at the first and, respectively, 
the second factor in the integrand in (fj) as follows: 



k 2 - 2 Pi k 



E 

j=0 



2V 



1,2. 



(4) 



All the resulting one- loop integrals are easily evaluated for any order of the ex- 
pansion. The operator A^o gives the following contribution at e ^ 0: 



M F r = -iTT 



d/2_ 



1+6 



E 

n=0 



r(n + 1 + e)r(-n - e) 
r(l -n -2e) 



g 2 



(5) 



The terms Ali-Fr and ftA 2 Fx are not individually regulated by dimensional regular- 
ization but their sum exists for general e ^ 0: 



(A4i + M 2 )F T = in 



d/2_ 



q 2 (m 2 ) e 



r(e)r(i-e)X: 



n=0 



n!r(l — n — e) 



\n(Q 2 /m 2 ) + if>(e) + ip{n + 1) - - e) - - n - e) 



m 



le=0 



<7 — > — OG 



27T 

g 2 



1 



1 



Li 2 - + - hr t - In t ln(t - 1) - -vr 



where Li 2 is the dilogarithm and t = Q 2 /m 2 . 



(6) 



Here ip is the logarithmic derivative of the gamma function. In the limit e — ► 0, one 
arrives at the known result: 



(7) 



3 Diagram with one non-zero mass 

The two-loop Feynman integrals under consideration can be written as 

F i (p 1 ,p 2 ,m,e) 

d d kd d l 

(I 2 - 2 Pl l)(l 2 - 2p 2 l)(k 2 - 2 Pl k)(k 2 - 2p 2 k){k 2 - m 2 5i )((k - I) 2 - m 2 ) ' ^ 

where m 51 = and m 52 = m. Let us apply the above prescriptions to them. The 
family of subgraphs that give non-zero contributions is shown in Fig. 2(a-f). Let 
us begin with the first case, when subgraphs (c) and (d) also happen to give zero 
contributions. 

The contribution of Fig. 2a is obtained by expanding the massive propagator in 
the geometric series in m. In the leading order, this is Fig. 2a with pure zero masses 
which was calculated in [[10], [□], [12| : 

(m d l 2 \ 2 l r 1 

[-G 2 (2,2)G 3 (2 + e,l,l)-G 2 (2,l) 



(Q 2 ) 2+2e e 



^G 3 (2,l,l + e) + G 3 (l,l,l) 



(9) 
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where G2{a,b) and Gs(a, b, c) are expressed through gamma functions and give, re- 
spectively, the value of the general one-loop propagator-type integral and triangle 
integral with the given kinematics, in the case of arbitrary indices of lines. 

All the resulting Feynman integrals in arbitrary order of the expansion are evalu- 
ated by integration by parts [13j (as in [|T^j) and are expressed, for general e, in terms 
of gamma functions. 

To write down the contribution of Fig. 2b it is necessary to expand the initial 
diagram, with the 6th line omitted, in Taylor series (under the sigh of the integral 
over the loop momenta) in the loop momentum flowing through this line in the whole 
diagram and then insert the result into the reduced diagram which is the tadpole 
generated by the 6th line. This subgraph Fig. 2b contributes to the expansion start- 
ing from the next-to-leading order, m 2 /(Q 2 ) 3 . The calculation of this contribution 
is rather simple and reduces to successive application of the one-loop integration 
formulae, with the G-functions G 2 and G 3 in the result. 

What is more non-trivial is the evaluation of the contributions of the subgraphs of 
Fig. 2e and f. The contribution of Fig. 2e is obtained by expansion of the propagators 
l/(/ 2 — 2 Pl l) and 1/ (k 2 — 2pik) in geometrical series in I 2 and k 2 , respectively: 

v rr d d kd d i(~ir^(i 2 )^ 

jfaj J (-2 Pl iy+ N i(l 2 -2p 2 l)(-2 Pl ky+ N 3(k 2 -2p 2 k)(k 2 y- N 3((k-l) 2 -m 2 )' [ 1 

However these contributions are not individually regularized by dimensional regular- 
ization so that let us introduce, temporarily, analytic regularization into the lines 
number 3 and 4, with A 3 7^ A 4 , then calculate the integrals involved and switch off 
the regularization in the sum. 

Thus the problem reduces to calculation of the following family of integrals: 



J(N 1: ...,N 7 



d d kd d l{l 2 ) N ' 



-2 Pl iy+ N ^(i 2 - 2 P2 iy+ N 2(-2 Pl ky+ N 3(k 2 - 2 P2 ky+ N * 

1 



(k 2 -m\) N s{{k-l) 2 -m 2 )^ 



(11) 



with integer N 6 and N 7 . Using integration by parts |l3j it is possible to express any 



integral with N 7 > through a linear combination of integrals with N 7 = 0. 

For integral (|II|) with N-? = 0, let us use the a-parametric representation. Then, in 
the resulting 6-fold integral, it is possible to apply the Mellin-Barnes representation in 
such a way that the internal integration over a-parameters can be performed explicitly 
and we are left with the following representation: 

(m d / 2 ) 2 

J(N u N 2 ,N 3 + \ 3 ,N 4 + \ 4: N 5 ,N 6 ,0) 



^Q2y+N 1 +N 3 +X 3 ^ m 2ye+N 2 +N 4 +N 5 +N( i +\ 4 

r(-e -N 3 -N 5 - A 3 )r(-e - N 4 - N 5 - \ 4 )T(2e + N 2 + N 4 + N 5 + N 6 + A 4 ) 

r(i + jv 2 )r(i + n 4 + A 4 )r(i + n 5 )t(i - e - Ar 5 )r(-Ar 1 - iv 3 - e - a 3 ) 
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F(s - e - N 2 + I) 



xr(s - N x - N 3 - e 



'r( s -e-JV3-A 3 + l) 
A 3 )r(e + N 2 + N 4 + A 4 - s)r(-s) . (12) 



The contour of integration is chosen in the standard way: the IR poles are to the 
left and the UV poles are to the right. The singularities at A3 — > A4 are due to the 
fact that some UV poles become identical with IR poles. By shifting the contour of 
integration these pole contributions are explicitly picked up. The rest of the integral 
is finite in the limit A 3 , A 4 — > 0, and the resulting integrals are easily evaluated by 
recurrence relations and finally using well-known integrals with four gamma functions. 

Following this procedure, any integral that is present in the contribution of Fig. 2e 
and f can be analytically evaluated, for arbitrary e. For example, the leading contri- 
bution takes the form 



U7T 



d/2\2 



(Q 2 ) 2 (m 2 ) 2e 



r(-e)r(e)r(2e)[L + ^(2e) + ^(1 + e) - 2^(1)] 



(13) 



where L = lnQ 2 /m 2 . To avoid the Euler constant in the results we present several 
first terms of the expansion of the function F 1 ■ (Q 2 ) 2 / (m d/2 Y(l + e)) 2 = J2n=o C n t~ n : 



Co 
Ci 



2 L +3^ 



1 1 T 4 1 1 ^2 t 2 

-(L + l)j,+ 

2 T 3 



I r3 
3 L 



WL - 3C(3) 



1 K(3)/- • rr- 



61 
360 

L 2 -L + ln 2 ^ 1 



-|L d + 3L Z - 6L 



■L - 26 + W + 11C(3) , 



= -&+i)$ + (lL*-lL-i + ± 



7T" 



_i L 3 + m L 2 _ 5L 



7T 2 L 
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8 



12 
17 



^3 = -(l^+l)i + (^ 2 -i^-! + > 



2 y + ¥c(3) 



-§L3 + iil 2 - f l - y L - m + ¥ y + ^ c(3 ) . 



(14) 



These results have been obtained with the help of the Mathematica system [T3| . 



4 Diagram with two non-zero masses 

In the case of the second Feynman integral that we consider there are contributions 
from all the subgraphs of Fig. 2. The contribution of Fig. 2a is obtained by expanding 
the propagators of the 5th and the 6th lines in the geometric series in m. As in the 
previous case the resulting integrals are calculated by recurrence relations based on 
integration by parts. The leading contribution has the same form The contri- 
bution of Fig. 2b in arbitrary order also poses no difficulties, with a result expressed 
through the G-functions G 2 and G 3 . 
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Consider now Fig. 2c and d. According to the general prescription, the contri- 
bution of Fig. 2c is written through the product of the expansion of the propagator 
l/(/c 2 — 2pik) in geometrical series with respect to k 2 and an expansion of the lower 
triangle subdiagram which is a function f(qf, q%, qf) of the following external momenta 
squared: q\ = k 2 — 2p±k, q\ = k 2 — 2p 2 k and q 2 = q 2 . So, this second expansion of / 
is in Taylor series (again in the sense of the expansion under the sign of integral) in 
g|, with the subsequent expansion of the result, as a function of k 2 — 2p\k in k 2 . All 
the derivatives of that triangle in q\ at q\ = are calculated by recurrence relations 



based on integration by parts (see, e.g., results for the first two derivatives in |14|). 

This product of the expansion of the 4th line and the lower triangle is then inte- 
grated with the product of other propagators (number 4 and 6), using simple one-loop 
integration formulae. However this contribution of Fig. 2c alone is not regularized 
dimensionally so that it is natural to consider the sum of both contributions of Fig. 2c 
and d. Still to handle these terms individually, one can introduce a temporal ana- 
lytic regularization, into the lines number 3 and 4, calculate them and switch off the 
analytic regularization in the sum. 

Following this procedure, the contribution of Fig. 2c and d can be evaluated 
in terms of gamma functions in every order of the expansion, for arbitrary e. For 
example, the leading, 1/(Q 2 ) 2 , contribution takes the form 

(iix d l 2 ) 2 r(e) 2 r(l -ef f Y(-ef , , . ,A 

' 1 -2 w l o \ +lnt + ^(e)-2^(-6) + ^(l) . (15) 



(Q 2 ) 2 +^m 2 Y eT(l-2e) [ T(-2e) 

To calculate the contribution of Fig. 2e and f it is reasonable to use the Mellin- 
Barnes representation for the 5th propagator and reduce the problem to calculation 
of the corresponding contribution for the first of our integrals, in the case when this 
line is analytically regularized. In particular, the leading, 1/(Q 2 ) 2 , contribution is 
obtained as the following integral of the analytically regularized version of eq. ([13]): 



IZ7T 



d/2\2 i 



J dsr(2e + s)T(e + s)T(-e - s)T(-s) 



(Q 2 ) 2 +*(m 2 y 2m Jc 

x [hit + i)(2e + s) + 4){l + e + s) - 2^(1 + s)) . (16) 

The contour C is chosen in the standard way, with the qualification that the pole 
s = — e (which is both UV and IR) is to the right of it. The pole and the finite part 
of eq. (|16D are evaluated separately and give the following result written up to e°: 

Hn d / 2 ) 2 f_, . . ^ 



{2C(3)lnt+^ 



(Q 2 ) 2 (m 2 ) 2e I s w 30 
+r(-e)r(e)r(2e) [hit + V(l - e) - 2^(1 - 2e) - ip(-e) + ip(e) + ^(2e)]} . (17) 

Now, the leading order of the asymptotic expansion of the second of our Feyn- 
man integrals is given explicitly by the sum of three terms corresponding to Fig. 2a, 
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Fig. 2(c&d) and Fig. 2(e&f), respectively, (g), QT5|) and ([17]). Note that the first 
term involves infrared and collinear divergences, the third term involves ultraviolet 
and collinear divergences, and the second term possesses all the three kinds of the 
divergences. However all the poles, up to 1/e 4 , are canceled in the sum and we get 

F 2 (Q 2 ,m 2 ) Q2 -°° 7^f^ln 4 t+^ln 2 t-6C(3)lnt+ 3l7r '^ 



(Q 2 ) 2 \2A 3 * w 180 J 

These results are in a good agreement with numerical MC calculations: e.g., there 
is 2% accuracy at t — 50. However this accuracy is achieved if all the terms in the 
parenthesis are taken into account. In other words, it is necessary to exhaust all 
the powers of the logarithm until the remainder is determined by the next power, 
m 2 /(Q 2 ) 3 . 

The calculation of the contribution of Fig. 2e and f in arbitrary order of the 
expansion happens to be more elaborated and will be reported in future publications 
as well as details of the present calculations, similar results for non-planar diagrams, 
with subsequent application of the new information obtained in two loops to extend 
well-known results on asymptotic behaviour in the Sudakov limit in QED and QCD 
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Figure 1: (a) One- loop vertex diagram, (a)-(c) Subgraphs contributing to the asymp- 
totic expansion of the diagram (a) in the Sudakov limit. 




(e) (/) 

Figure 2: (a) A typical two-loop vertex diagram, (a)-(f) Subgraphs contributing to 
the asymptotic expansion of the diagram (a) in the Sudakov limit. 
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